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We consider a generalization of the contact process stochastic model, including an additional 
autocatalitic process. The phase diagram of this model in the proper two-parameter space displays 
a line of transitions between an active and an absorbing phase which starts at the critical point of the 
contact process and ends at the transition point of the voter model. Thus, a crossover between the 
directed percolation and the compact percolation universality classes is observed at this latter point. 
We study this crossover by a variety of techniques. Using supercritical series expansions analyzed 
with partial differential approximants, we obtain precise estimates of the crossover behavior of the 
model. In particular, we find an estimate for the crossover exponent 4> — 2-00 ± 0.02. We also show 
arguments that support the conjecture (f> = 2. 
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I. INTRODUCTION 



The phase transitions exhibited by stochastic models with absorbing states have attracted much attention in recent 
years, particularly in order to identify and understand the aspects which determine the universality classes in those 
models. Most of these models have not been solved exactly, but a variety of approximations allow quite conclusive 
results regarding their critical properties. Stochastic models are, of course, well fitted for simulations, but closed form 
approximations and other analytical approaches have also been useful in investigating their behavior 

One of the simplest and most studied model of this type is the contact process (CP), which was conceived as a 
simple model for the spreading of an epidemic and proven to display a continuous transition between the absorbing 
and an active state, even in one dimension Q. Actually, it was found that the CP is equivalent to other models 
such as SchlogPs lattice model for autocatalytic chemical reactions Q and Reggeon Field Theory (RFT) Q. The CP 
belongs to the direct percolation (DP) universality class, together with others models such as the Ziff-Gulari-Barshad 
model of catalysis Q and branching and annihilating walks with an odd offspring [6]. The DP conjecture states that 
all phase transitions between an active and an absorbing state in models with a scalar order parameter, short range 
interactions and no conservation laws belong to this class 0|. This conjecture was verified in all cases studied so far 

@ 

Here we study a generalization of the CP, with an additional parameter, so that the CP transition point becomes a 
critical line. Since the symmetry properties of this generalized model are the same of the CP, it is expected that this 
critical line should belong to the DP universality class. However, at one point of this line the model is equivalent to 
the zero temperature Glauber model [jj, also called the voter model [Tll j. which displays a spin inversion (or particle- 
hole) symmetry and therefore belongs to another universality class. Thus the critical line in the phase diagram of 
the generalized model starts at the CP model and ends at the voter model, a crossover between the two universality 
classes being observed. The voter model belongs to the compact percolation universality class, and also corresponds 
to a limiting point in the phase diagram of the Domany-Kinzel cellular automaton, where an exact solution is possible 
Thus, the exact critical exponents are known for this model. In a study of models with several absorbing states 
simulational results are shown for the model we are considering here and a study of the crossover between direct 
percolation and compact direct percolation may be found in 0], motivated by the possibility to explain the non- 
universality in models with several absorbing states as a surface effect. AlsOjthe shape of the critical line close to the 
CDP endpoint in the Domany-Kinzel automaton was studied in detail [IH , \wi \vfy . and these results are compared to 
our findings in the conclusion. Some physical motivation for the model we are studying here might be given. In the 
contact process, the additional term might be understood as an enhancement of the possibility of a sick individual 
to recover proportional to the number of its first neighbor which are healthy. However, our motivation to study the 
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model is centered on its simplicity and the universality class crossover present in its phase diagram. 

In section [H] we define the model and explain how supercritical series expansions may be obtained for it. The 
coefficients of the two- variable series for the survival probability up to order 25 are given. Section IIIII contains the 
description and the results of the Pade and PDA estimates for the model, with emphasis on the multicritical behavior 
in the voter model limit. Final discussions and the conclusion may be found in section llVl 



II. DEFINITION OF THE MODEL AND CALCULATION OF THE COEFFICIENTS OF THE 
SUPERCRITICAL SERIES FOR THE SURVIVAL PROBABILITY 

The model is defined on a one-dimensional lattice with N sites and periodic boundary conditions. Each site is 
occupied either by a particle A or a particle B, no holes are allowed. The microscopic state of the model may thus be 
described by the set of binary variables r\ = (r/i, 772, ■ ■ ■ , J?jv), where r/i = or 1 if site i is occupied by particles B or 
A, respectively. 

The model evolves in time according to the following Markovian rules: 

1. A site i of the lattice is chosen at random. 

2. If the site is occupied by a particle B, it becomes occupied by a particle A with a transition rate equal to p a n,A/2, 
where ua is the number of A particles in the sites which are first neighbors to site i. 

3. If site i is occupied by a particle A, it may become occupied by a particle B through two processes: 

• Spontaneously, with a transition rate p c . 

• Through an autocatalytic reaction, with a rate p b riB/2, where ub is the number of B particles in the sites 
which are first neighbors to site i. 

We define the time in such a way the the non-negative parameters p a ,Pb, and p c obey the normalization p a +Pb +Pc = 
1. We may then discuss the behavior of the model in the (p a ,Pc) plane without loss of generality. 
The probability P(rj,t) to find the system in state 77 at time t obeys the master equation 

= EK(V)P(V,t) - wMP(v,t)} (i) 

i=l 

where rf corresponds to the following configuration 

= (171,..., l-rfe,...,7?jv) (2) 
and Wi(j]) is the transition rate of the model, given by 

Wi(v) = ^^-JVij^Vi+s + Vi, (3) 

1 8 

where fj, = p a /(l — p a ), 7 = (1 — Pc)/Pa, and the sum is over first neighbors of site i. 

It may be useful to remark that this model may be mapped to a spin system if we describe sites occupied by A and 
B particles by an Ising spin variables <7, = 1 and Oj = —1, respectively. In these variables, the transition rate will be 
given by 
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(4) 



where a = (p a +p b + 2p c )/2, (3 = (p a -p b -2p c )/(pa+p b + 2p c ),e = (p a + p b ) / \p a + p b + 2p c ) , and £ = (p a -pb)/(j>a + 
p b + 2p c ). 

In two particular cases, this model corresponds to well known models. If we make p b = or 7 = 1 contact process 
is recovered Q: 

) ('7) = f(l-^)E^ + 7 ?- ( 5 ) 



If now we take p a = Pb and p c = in the spin formulation of the model, the zero temperature linear Glauber model 
[Icj . also known as the voter model, is recovered [Til ]: 



(LGM), \ a 
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(6) 



This model has been studied using mean-field approximations ^3], as we ll as simulations fili. Hflj . For p c > 
0, the stationary state at low values of pa corresponds to the absorbing state, where the density of A particles 
Pa =< Na > /N vanishes. As p a is increased, a continuous phase transition occurs and an active stationary 
state (pa > 0) is stable at high values of y„,. Thus a critical line is present in the phase diagram starting at 
(p a = 0.767325(6), p c =l-p a = 0.232674(4)) [1.9j (contact process), and ending at (p a = 1/2, p c = 0) (linear Glauber 
model), where the transition is discontinuous and between two absorbing states (pa = 1 and pb = 1). In figure 
results from mean-field calculations and simulations for the phase diagram are shown 0|. While it is expected that 
the critical exponents at the whole critical line are the ones of the directed percolation (DP) universality class, a 
crossover to the compact directed percolation (CDP) universality class exponents happens as p c vanishes yj. Thus, 
we may recognize the point (p a = 1/2, p c = 0) of the phase diagram as a multicritical point. Using supercritical series 
expansions, for the survival probability, we will study the multicritical singularity in the neighborhood of this point. 

Now let us develop a two-variable supercritical series expansion for the model. We follow closely the operator 
formalism presented in the paper by Jensen and Dickman on series for the CP process and related models We 
may then represent the microscopic configurations of the lattice by the direct product of kets 

\V) = <S> \Vih (7) 

i 

which are defined to be orthonormal 

(vW)=hk,w (8) 

i 

Now we may define A particles creation and annihilation operators for the site i: 

A\Vi) = (l-»fc)|»Ji + l), 

Mm) = m\m - 1)- (9) 

In this formalism, the state of the system at time t is 

m)) =1>(»7,*)|»7>. (10) 
M 



If we define the projection onto all possible states as 



the normalization of the state of the system may be expressed as ( \ip) — 1. In this notation, the master equation for 
the evolution of the state of the system (Eq. is 

= s\m). (12) 

The evolution operator S may be expressed in terms of the creation and annihilation operators as S — XSq + V where 

S a = ^[a(2-4_ 1 A i _ 1 +4 +lJ 4 i+1 ) + l](A i -4A i ), (13) 

i 

V = £(4 + A\A t - l)(4-i^-i + A\ +1 A i+1 ), (14) 



and the new parameters 



A = — 

Pa 



and 



Pb 

a = - — 

2p c 

were introduced. 

We notice that the operator So only annihilates A particles (transitions A — y B), while the operator V acts in the 
opposite way, generating transitions B — > A. Thus, for small values of the parameter A the creation of A particles is 
favored, and the decomposition above is convenient for a supercritical perturbation expansion. Let us show explicitly 
the effect of each operator on a configuration (C). 



S (C) = aJ2(C'i) + 2a£(Cj) + E(^) - + + r K C ), 



(15) 



where the first sum is over the r\ sites with A particles and one B neighbor, the second sum is over the r 2 sites 
with A particles and two B neighbors and the third sum is over all r sites with A particles of the configuration (C). 
Configuration (C'A is obtained replacing the A particle at site i by a B particle. The action of operator V is 



V(C) = + 2$>") - (qi + 2q 2 )(C) 1 



(16) 



where the first sum is over the q\ sites with B particles and one A neighbor, the second sum is over the q 2 sites with 
B particles and two A neighbors. Configuration (Cf) is obtained replacing the B particle at site i in configuration (C) 
by a A particle. 

To obtain a supercritical expansion for the ultimate survival probability of A particles, we start by remembering 
that in order to access the long-time behavior of a quantity, it is useful to consider its Laplace transform. For instance, 
the Laplace transform of the state of the system is 



1^)) = / e- st m)), 

Jo 

and inserting the formal solution \ip{t)) = e st |-0(O)) of the master equation El we find 

$( S )) = (s-s)- i mo)}- 

The stationary state |^(oo)) = lim^oo \ip(t)) may then be found noticing that 

ji/K 00 )) = lim s|-0(s)), 



(17) 



(18) 



(19) 



which may be obtained integratingElby parts. A perturbative expansion may be obtained assuming that \tp(s)) may 
be expanded in powers of A and usinglTsl 
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(20) 



(21) 



(22) 



(23) 



The action of the operator (s — V) 1 on an arbitrary configuration (C) may be found noting that 



(s-V)- 1 (C) = s~ 1 (C) + 



V 



s(s - V) 

and using the expression ITfil for the action of the operator V, we get 



( S -V)- 1 (C) = s q {(C) + ( S -V)- 1 



XX) + 2£(C?) 



(24) 



(25) 



where the first sum is over the q\ sites with B particles and one A neighbor, the second sum is over the qi sites with 
B particles and two A neighbors, and we define s q = l/(s + q), where q = qi + 2q 2 . 

It is convenient to adopt as the initial configuration a translational invariant one with a single A particle (periodic 
boundary conditions are adopted). Now we may notice in the recursive expression |2£l that the operator (s — V)^ 1 
acting on any configuration generates an infinite set of configurations, and thus we are unable to calculate in a 
closed form. We may, however, calculate the extinction probability p(s), which corresponds to the coefficient of the 
vacuum state |0). As happens also for the models related to the CP studied in ^| configurations with more than j 
particles only contribute at orders higher than j, and since we are interested in the ultimate survival probability for 
A particles = 1 — lim s _,o sp(s), s q may be replaced by 1/q in Eq. |23 To illustrate the procedure, we will perform 
the explicit calculation of the series for lim s ^ sp(s) up to third order in A. We furthermore represent a configuration 
denoting by • a site occupied by an A particle and by o a site occupied by a B particle. Thus (• o •) denotes the 
translationally invariant configuration J^i A\a\ +2 1 0). The B particles situated to the left of the leftmore A particle in 
a configuration and the ones situated to the right of the rightmost A particle are omitted in the representation of this 
configuration. The vacuum state will be represented by (0). As stated above, at t = the system will be supposed to 
be in the configuration \ip{0)) = X^^ll^) = (•)■ Keeping configurations with up to three particles, and omitting the 
global factor 1/s, the first of the recursion relations Eqs. [22] leads to 



l^o) = 5 [(•) + (") + (•«•) + 



(26) 



The next step is the calculation of \<po) — S'olV'o}- Now we need to keep only configurations with up to two particle. 
The result is 



|0o) = i[(2a+l)(O) + (.) + .o.) + 
Now we obtain l^i) = (s — V)^ 1 ]^) for s = 0, resulting in 

From this point on we will only show the results of each step, up to the relevant numbers of A particles: 

\4>i) 



i(2a + l)(0) + i(.) + i(..) + i(.o.) + 



(27) 



(28) 



I(2a+l)(0) + i(2a + 2)(«) + 



\i> 2 ) = i(2a+l)(0) + i(2a + 2)(.) + " 
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The first coefficients of the ultimate survival probability will then be given by 

Poo = 1 - l(2a + 1)A + --(2a + 1)A 2 + \{2a + 2)(2a + 1)A 3 + 
2 4 o 



(29) 



The algebraic operations above may be easily performed in a computer using a proper algorithm. The configurations 
are expressed as binary numbers and the coefficients as double precision variables. Although we have tried to do the 
calculation representing the coefficients as rational numbers, thus avoiding any roundoff errors since all calculations 



were done with integers, we found that the denominators increase very rapidly with the order of the calculations, and 
thus we were unable to perform the calculations this way up to a reasonable order. With rather modest computational 
resources (Athlon MP2200, double processor, 1Gb memory) it is not difficult to calculate the coefficients up to order 
25. The required processing time amounts to about 6 hours, the limiting factor is actually the memory required for 
the calculation. The maximum number of terms (polynomials in a) amounts to more than 44 x 10 6 . We define the 
coefficients bij as: 

1 25 i-1 

Poo = 1 - -(2a + 1)A Y.Y. h ' ' X ' (r '- ( 3 °) 

and they are listed in Table |U Up to order 24, our results are numerically coincident with the supercritical series 
expansion for the ultimate survival probability of the contact process ^| , in the particular case a = (we remark 
that the variable A in the supercritical expansion for the ultimate survival probability in reference 0] is half the 
variable A we use here). 



III. ANALYSIS OF THE SERIES 



Let us consider initially the one variable series for fixed values of a 

25 

Poo = ^ ai(a)X i . (31) 

i=0 

As a preliminary test, we may apply the ratio method j^l to these series. The results for = ai/a.;_i as functions 
of 1/i are depicted in FigureEJ In the case a = (circles), the asymptotic linear behavior r.; « 1, 6489(1 — 0.7231/i), 
obtained from precise estimates for A c and the exponent (3 for the contact process (Tg j. In the figure it is apparent 
that the ratios approach the asymptotic limit as i is increased, as a matter of fact this approach is close to linear 
in 1/i. Thus, we may infer that for this model the singularity which is closer to the origin is actually the physical 
singularity. As the value of a is increased, one may notice that the asymptotic linear behavior in 1/i for the ratios 
occurs only for higher values of i, and thus it will be increasingly difficult to obtain precise estimates for higher values 
of a. 

Even in cases where the singularity of physical interest is the one closest to the origin, the d-log Pade approximants 
usually lead to better estimates than the ratio method [20] . The approximants are defined as ratios of two polynomials 
P L (X) and Q M (A): 

Flm{X) ~ cma) ~ TT^F (32) 

The series for ^lnP oc (A) for fixed values of a are substituted in the defining equation and the coefficients of 
the polynomials are chosen such that the identity is true up to the order of the available series expansion. Thus 
approximants with L + M < 24 may be built with the available series. Usually diagonal (L = M) or close to diagonal 
approximants furnish better results, so we restricted our calculations to these cases. The estimate for the critical 
value of A is found among the poles of the approximant, the estimate for the critical exponent (3 will be the residue 
at this pole. 

We thus built approximants with a ranging between and 40, estimating the critical value of the parameter A 
as well as the critical exponent (3. Although for small values of a the results are very good, with estimates of 
comparable to the best ones in the literature for the contact process, as the value of p c is decreased we notice a 
growing dispersion of the estimates for the exponent and for p c < 0.01 most of the approximants lead even to ill 
conditioned systems of linear equations for the coefficients, and therefore we were not able to obtain estimates in this 
region. We made additional one- variable investigations, such as obtaining approximants for (A c — A)^P,5(A) = (3 for 
several values of A c and searching for the intercept of the curves (3(\ c ) and non-homogeneous Pade approximants 
[2fT.l2lj. and although some improvements of the estimates may be obtained in certain cases, the situation does not 
change qualitatively. In figure the Pade estimates for the critical line are displayed. 

Actually, the increasing dispersion of the estimates as the parameter p c approaches zero is not a surprise, since as was 
mentioned above in this limit the model corresponds to the voter model which is in the compact directed percolation 
(CDP) universality class, whose exponents are different from the ones of the contact process, which belongs to the 
directed percolation (DP) universality class. For the voter model, the exponent (3 of the order parameter is equal to 



zero (a discontinuity in the order parameter occurs at the transition) pj, but the exponent for the survival probability 
j3' is distinct from (3 and equal to 1 [13] ■ Therefore a crossover to from the DP to the CDP universality class occurs 
as p c — > 0, and it is known that in such situations the reduction of two- variable series to one variable leads to very 
poor estimates So we analyzed the series without reducing the problem to one variable, and to our knowledge 
the best results for two-variable series ap plie d to multicritical phenomena in the literature were obtained using the 
partial differential approximants (PDA's) [23]. They may be regarded as a generalization to two variables of the d-log 
Pade approximants. The defining equation of the approximants is 

PL(x, y)F(x, y) = Q M (x, v)^^ + v) ^ 2 ^ , (33) 

where P, Q, and R are polynomials in the variables x and y with the set of nonzero coefficients L, M, and N, 
respectively. The coefficients of the polynomials are obtained through substitution of the series expansion for the 
quantity which is going to be analyzed 

f(x,y)= V (34) 

k,k'=0 

into the defining equationESIand requiring the equality to hold for a set of indexes defined as K. Again this procedure 
leads to a system of linear equations for the coefficients of the polynomials, and since the coefficients fk,k' of the series 
are known for a finite set of indexes this sets an upper limit to the number of coefficients in the polynomials. Since the 
number of equations has to match the number of unknown coefficients, we must have that the numbers of elements 
in each set satisfy K = L + M + N — 1 (one coefficient is fixed arbitrarily). An additional issue, which is not present 
in the one-variable case, is the symmetry of the polynomials. Two frequently used options are the triangular and the 
rectangular arrays of coefficients. The choice of these symmetries is related to the symmetry of the series itself [24l |. 
Below we discuss the solution we adopted in the present case for this point. 

Let us suppose that quantity represented by the series is expected to have a multicritical behavior at a point (x c , y c ), 
described by 

f(x,y)^\Axrz(^}-), (35) 



where 

Ax = (x-x c )-(y-y c )/e2, (36) 

and 

Ay = (y - Vc) - ei(x - x c ). (37) 

Here 7 is the critical exponent of the quantity described by / when Ay = 0, e\ and e-2. are the scaling slopes and 
4> is the crossover exponent. The function Z(z) is singular for one or more values of its argument, corresponding to 
the critical line(s) incident on the multicritical point. Once the coefficients of the defining polynomials are obtained, 
the estimated location of the multicritical point corresponds to the common zero of the polynomials Qm and Rn- 
This may be seen substituting the scaling form EH in the defining equation ESI of the approximant. The exponents 
and scaling slopes may also be obtained directly from the polynomials, without integrating the partial differential 
equation. A detailed discussion of the algorithm, as well as computer codes, may be found in j24[. 

Before proceeding with the analysis of the series, it is convenient to perform a change of variables, since the 
multicritical point in the original variables is located at a — * 00. We thus express the series in the variables 

x = \= 2 -^ 



y = a\ 



Pa 

1-Pa-Pc 
Pa 



In these new variables the multicritical point is located at x = 0, y = 1, and the survival probability may be written 

as 



Poo = 1 - t, x - y ~ xF(x,y), 



(38) 



where 



24 i 

Fir.,,; XX'''- 1 -'-' V ( 39 ) 

i=l j=0 

is represented by a series with triangular symmetry, which is very convenient to be analyzed using PDA's. The 
number of approximants which may be obtained from the series is very big, so we restricted ourselves to approximants 
with the number of elements in M close to the number of elements in N. The polynomials had the same triangular 
symmetry as the series, but in most cases some higher order elements of the polynomials were set equal to zero in 
order to match the number of unknowns to the number of equations in the set of linear equations for the coefficients. 
This is a rather standard procedure and is discussed in detail by Styer 0|. Even at rather low orders, we found a 
reasonable agreement between most of the estimates from different approximants. Finally, we considered a set of 42 
approximants which use the elements of the series for F(x, y) with highest orders i between 15 and 24. Discarding 
some approximants which generated estimates which were rather far away of the general trend, finally we used a set 
of 36 approximants to obtain the estimates. 

In figure 21 the estimates for the location of the multicritical point are displayed. We notice that the estimates are 
very close to the exactly known values x — and y = 1. The estimated values for x c and y c are (0.4 ± 1.8) x 10~ 6 and 
y = 1.0000 ±0.0003. The exponent —7 in the multicritical scaling form |2H corresponds to the exponent f3' of the CDP 
universality class. The estimated value is equal to 1.00 ± 0.01, which agrees very well with the exact value 0' = 1 
0. Finally, the crossover exponent was estimated as (j> = 2.01 ± 0.03, and thus the mean field value for this exponent 
{(f) — 2 [1^1) is within the confidence interval of our estimate. The estimates for /?' and <j) are shown in figure We 
also obtained biased PDA's, fixing the other parameters at their known values and calculating improved estimates 
for (j). This procedure resulted in the estimate (f> = 2.00 ± 0.02. We thus conclude that our estimates for cf) are very 
close to the classical value of the crossover exponent. The estimates for the slopes of the scaling axes show a rather 
broad distribution, a significant majority of the approximants provide quite large values for e±, while 62 is typically 
much smaller. This suggests that Ay = x and it is reasonable to choose Ax = l — y, since it corresponds to the weak 
direction, parallel to the critical line at the multicritical point. In the limit of the voter model (x — 0) we notice that 
the series EHI reduces to the exact result Poo = 1 — y. 

One way to actually estimate values of the quantity which is described by a PDA is to integrate the equation using 
the method of characteristics. A timelike variable r is defined and a family of curves (x(t), y(r)) (the characteristics) 
is considered. These curves are defined by the equations 

^ = Qm(x(t),v(t)), 

^ = R N (x(r),y(T)). (40) 

Along such a curve, the defining equation of the PDA EH leads to an ordinary differential equation for F. 

^ = P L ( x (r),y(r))F, (41) 

which may readily be integrated, providing the value of F at the points of the characteristics, once we know this value 
at a initial point. Our efforts to obtain estimates for the survival probability, particularly close to the multicritical 
point, using this procedure were not very successful. Nevertheless, it is worth to mention that an estimate for the 
critical line may be obtained this way. The critical line, which connects the point which corresponds to the CP to the 
multicritical point of the voter model transition is a line of singularities, and it is not difficult to show that such a line 
is one of the characteristics defined by equations Therefore, the characteristic whose initial point is located on 
the transition point of the CP, which is known with great accuracy, is an estimate for the critical line. We therefore 
chose a particular approximant with estimates close to the mean values. The number of elements in the sets for this 
approximant are K = 300, L = 171, and M = N = 65. In figure El this characteristic curve is shown, and a nice 
agreement with other estimates of the critical line may be observed. 



IV. CONCLUSION AND DISCUSSION 



We may notice that the coefficients b(i,i — 1), for i = 1,2,... are equal to 1/2, a result which is indicated by the 
coefficients below but may be shown to be true in general using the operator formalism above. Therefore, we may 



sum these set of terms in the series, obtaining 

Poo = (l - y) 



2(l-y)* 



(42) 



Now if we compare this expression with the multicritical scaling form 1351 we may recognize between braces the two 
first terms of a Taylor expansion of the multicritical scaling function Z(z), where the variable z is identified as rrz^Tp ■ 
This agrees with the estimates obtained from the PDA's. Moreover, remembering that for y = the series represents 
the CP, the function Z(z) may be recognized as the survival probability of the CP, which was studied in great detail 
flfl| and found to have a singularity at z = 0.6064 with the exponent f3' = j3 = 0.276486. As expected, the critical line 
is characterized by the same exponent of the CP, that is, is in the DP universality class. The critical line corresponds 
to z — zq, that is 

x = z Q (l-y) 2 , (43) 

and this curve is shown in figure 03 It is interesting to notice that the agreement of this estimate of the critical line 
with the other two, obtained from Pade and partial differential approximants, is quite good even far away from the 
multicritical point. The multicritical scaling form with the identification of the scaling variable z above is exact in 
the limit p c = (biased voter model) and reproduces the supercritical series expansion for the CP, when y = 0. It 
does, however, not reproduce the two-variable supercritical series expansion for the full model we considered here 
away from these limiting cases. 

The Domany Kinzel probabilistic cellular automaton [p2j m P af t of its two-dimensional phase diagram corresponds 
to a synchronous update version of the contact process, and as stated above has the synchronous update voter model 
as a endpoint of the critical line. It is believed, although to our knowledge not proven, that if in a particular model the 
update procedure is changed from synchronous to asynchronous, the critical exponents do not change. Some bounds 
for the critical line are presented in Although an upper bound for this line due to Liggett [2^| is quadratic close 
to the CDP endpoint, the lower bounds are linear, and thus the asymptotic behavior of the critical line is not fixed 
by those bounds. The critical line is studied in more detail by simulations and series expansions in ^j, and based on 
this results the authors conjectured a quadratic asymptotic behavior of the critical line, consistent with <f> = 2. Finally, 
a more detailed series analysis is done in [17|, but since one- variable Pade approximants were used, the results are not 
precise in the region close to the CDP point. Thus, there are indications that the crossover exponent has the same 
value in both models, and if these indications are correct, the invariance of this multicritical exponent with respect 
to the update procedure is verified in this particular case. 

In conclusion, the analysis of the series for the ultimate survival probability using PDA's lead to quite precise 
estimates for the multicritical exponents, and these estimates, together with the possibility to sum the terms of the 
two- variable series which are linear in x allowed us to conjecture the exact form of the multicritical scaling expression. 
The multicritical scaling function Z(z) is known as a series expansion up to order 25. 
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25 


23 





75000019616897146690D- 


-02 


25 


22 





34503863852011384949D- 


-04 


25 


21 





75555516468201169288D- 


-05 


25 


20 





96061396122189321999D- 


-06 


25 


19 





79292390146647733217D- 


-07 


25 


18 





45593449242937929000D- 


-08 


25 


17 





19157385641492643557D- 


-09 


25 


16 





60860240979789494986D- 


-09 


25 


15 





14987944522389162749D- 


-10 


25 


14 





29148047590757961700D- 


-10 


25 


13 





45381694251898121450D- 


-10 


25 


12 





57125179306761664221D- 


-10 


25 


11 





58520654724976619576D- 


-10 


25 


10 





48967910880832077324D- 


-10 


25 


9 





33495648941633358042D- 


-10 


25 


8 





18689764057712466183D- 


hlO 


25 


7 





84597759965914463009D- 


-09 


25 


6 





30764609873623349756D- 


-09 


25 


5 





88524517163033653588D- 


-08 


25 


4 





19689926825686265843D- 


-08 


25 


3 





32638605561858318182D- 


-07 


25 


2 





37947216762104512000D- 


-06 


25 


1 





27603268426234146559D- 


-05 


25 








94519638900755200694D- 


-03 



Table I: Coefficients bij of the supercritical 



series expansion for the ultimate survival probability of A particles. 



■ ■ ■ One-site mean field 

— Two-sites mean-field 

— Simulations 




Figure 1: Phase diagram of the model obtained from mean-field approximations and simulations. The physical region is 
Pa + Pc < 1 and the line p a + p c = 1 corresponds to the contact process. In the region labeled RI the absorbing phase (all 
particles of type B) is stable, while in the region RII an active phase is stable, with a nonzero density of particles A. 
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Figure 2: Ratio of successive coefficients of the series expansion in A for fixed values of a. The results shown are for a = 
(circles), a = 1, a = 2, a = 3 a = 5, and a = 10 (triangles pointing down). The solid line shows the expected asymptotic 
behavior for a = (contact process). 




Figure 3: Estimates for the critical line. Curve (a) shows results from simulations (triangles) and Pade approximants (circles). 
Curve (b) displays the characteristic which starts at the CP (full line) and simulation results (triangles). In curve (c) the full 
line is again the characteristic and the values which follow from the multicritical scaling form are represented by crosses. 
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Figure 4: Estimates for the location of the multicritical point provided by the set of PDA's. 
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Figure 5: Estimates for the values of the exponents j3' and (j>. 



